REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. _ 

2.  REPORT  TYPE 

New  Reprint 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

611102 


5f.  WORK  UNIT  NUMBER 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 

14.  ABSTRACT 
see  attached 


15.  SUBJECT  TERMS 

systematic  errors,  acoustic  tomography,  temperature  and  velocity  fields 

17.  LIMITATION  OF  15.  NUMBER 

ABSTRACT  OF  PAGES 

UU 

Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Vladimir  Ostashev _ 

19b.  TELEPHONE  NUMBER 
303-497-3712 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

UU 

UU 

UU 

7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

University  of  Colorado  -  Boulder 
3100  Marine  Street  Room  479 
572  UCB 

Boulder,  CO _ 80303  -1058 _ 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 

1 1 .  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

57813-EV.7 


5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 


4.  TITLE  AND  SUBTITLE 

Sonic  Anemometer  as  a  Small  Acoustic  Tomography  Array 


6.  AUTHORS 

Sergey  N.  Vecherin,  Vladimir  E.  Ostashev,  Christopher  W.  Fairall,  D. 
Keith  Wilson,  Ludovic  Bariteau 


3.  DATES  COVERED  (From  -  To) 

5a.  CONTRACT  NUMBER 
W91  INF- 10- 1-04 15 


1.  REPORT  DATE  (DD-MM-YYYY) 


Report  Title 

Sonic  Anemometer  as  a  Small  Acoustic  Tomography  Array 


see  attached 


ABSTRACT 


REPORT  DOCUMENTATION  PAGE  (SF298) 
(Continuation  Sheet) 


Continuation  for  Block  13 


ARO  Report  Number  5781 3.7-EV 

Sonic  Anemometer  as  a  Small  Acoustic  Tomogr 


Block  13:  Supplementary  Note 

©2013  .  Published  in  Boundary-Layer  Meteorology,  Vol.  Ed.  0  (2013),  (Ed.  ).  DoD  Components  reserve  a  royalty-free, 
nonexclusive  and  irrevocable  right  to  reproduce,  publish,  or  otherwise  use  the  work  for  Federal  purposes,  and  to  authroize 
others  to  do  so  (DODGARS  §32.36).  The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and 
should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other 
documentation. 


Approved  for  public  release;  distribution  is  unlimited. 


Sonic  Anemometer  as  a  Small  Acoustic 


Sergey  N.  Vecherin,  Vladimir 
E.  Ostashev,  Christopher  W.  Fairall, 
D.  Keith  Wilson  &  Ludovic  Bariteau 


Boundary-Layer  Meteorology 

An  International  Journal  of  Physical, 
Chemical  and  Biological  Processes  in 
the  Atmospheric  Boundary  Layer 


BOUNDARY-L 

METEOROLOG 


ISSN  0006-8314 

Boundary-Layer  Meteorol 
DOI  1 0.1 007/si  0546-01 3-9843-9 


VOLUME  148  No.  2  August  2013 


An  International  Journal  of  Physical,  Chemical  and 
Biological  Processes  in  the  Atmospheric  Boundary  Layer 


Editor:  J.R.  Garratt 


,*0*0(4* 

■jf  % 

■  1 

4? 


4cj  Springer 


ISSN  0006-8314 


^  Springer 


Your  article  is  protected  by  copyright  and  all 
rights  are  held  exclusively  by  Springer  Science 
+Business  Media  Dordrecht.  This  e-offprint 
is  for  personal  use  only  and  shall  not  be  self- 
archived  in  electronic  repositories.  If  you  wish 
to  self-archive  your  article,  please  use  the 
accepted  manuscript  version  for  posting  on 
your  own  website.  You  may  further  deposit 
the  accepted  manuscript  version  in  any 
repository,  provided  it  is  only  made  publicly 
available  12  months  after  official  publication 
or  later  and  provided  acknowledgement  is 
given  to  the  original  source  of  publication 
and  a  link  is  inserted  to  the  published  article 
on  Springer's  website.  The  link  must  be 
accompanied  by  the  following  text:  "The  final 
publication  is  available  at  link.springer.com”. 


Author's  personal  copy 


Boundary-Layer  Meteorol 
DOI  10.1007/sl0546-013-9843-9 

ARTICLE 

Sonic  Anemometer  as  a  Small  Acoustic  Tomography 
Array 

Sergey  N.  Vecherin  •  Vladimir  E.  Ostashev  •  Christopher  W.  Fairall  • 
D.  Keith  Wilson  •  Ludovic  Bariteau 


Received:  6  December  2012  /  Accepted:  18  July  2013 
©  Springer  Science+Business  Media  Dordrecht  2013 


Abstract  The  spatial  resolution  of  a  sonic  anemometer  is  limited  by  the  distance  between  its 
transducers,  and  for  studies  of  small-scale  turbulence  and  theories  of  turbulence,  it  is  desirable 
to  increase  this  spatial  resolution.  We  here  consider  resolution  improvements  obtainable  by 
treating  the  sonic  anemometer  as  a  small  tomography  array,  with  application  of  appropriate 
inverse  algorithms  for  the  reconstruction  of  temperature  and  velocity.  A  particular  modifica¬ 
tion  of  the  sonic  anemometer  is  considered  when  the  number  of  its  transducers  is  doubled 
and  the  time- dependent  stochastic  inversion  algorithm  is  used  for  reconstruction.  Numerical 
simulations  of  the  sonic  anemometer  and  its  suggested  modification  are  implemented  with 
the  temperature  and  velocity  fields  modelled  as  discrete  eddies  moving  through  the  sonic’s 
volume.  The  tomographic  approach  is  shown  to  provide  better  reconstructions  of  the  tem¬ 
perature  and  velocity  fields,  with  spatial  resolution  increased  by  as  much  as  a  factor  of  ten. 
The  spatial  resolution  depends  on  the  inverse  algorithm  and  also  improves  by  increasing  the 
number  of  transducers. 

Keywords  Acoustic  tomography  •  Small-scale  turbulence  •  Sonic  anemometer 


1  Introduction 

Sonic  anemometers  are  robust  instruments  for  measurements  of  temperature  and  velocity 
and  have  been  used  for  such  measurements  for  many  decades,  e.g.,  Kaimal  and  Finnigan 
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(1994).  Due  to  concerns  about  flow  distortion,  the  transducers  of  a  sonic  anemometer  are 
located  at  some  distance  from  each  other,  and  as  a  result  a  sonic  anemometer  enables  only 
path-averaged  measurements  of  temperature  and  velocity.  Horst  and  Oncley  (2006)  reviewed 
currently  used  sonic  anemometers  and  showed  theoretically  that  such  an  instrument  enables 
the  correct  measurement  of  spectra  up  to  a  wavenumber  that  is  about  the  inverse  of  the  dis¬ 
tance  between  its  transducers.  There  are,  however,  several  important  applications/concerns 
in  boundary-layer  physics  and  turbulence  theory  that  require  analysis  of  turbulent  fields  at 
smaller  scales  (e.g.,  Fairall  et  al.  2003;  Chamecki  and  Dias  2004).  Among  these  are:  (i)  stud¬ 
ies  of  the  inertial  subrange  at  small  scales  that  are  important  for  analysis  of  the  dissipation 
rate  and  the  turbulent  kinetic  energy  budget,  (ii)  studies  of  turbulence  closure  models,  (iii) 
studies  of  energy  transfer  in  the  atmospheric  boundary  layer  that  are  important,  for  exam¬ 
ple,  for  wind  energy  applications,  (iv)  measurements  of  turbulence,  particularly  momen¬ 
tum  and  heat  fluxes,  when  the  energy-containing  range  extends  to  spatial  scales  smaller 
than  those  resolved  by  currently  used  sonic  anemometers,  e.g.,  near  the  surface  or  within 
canopies,  (v)  studies  of  small-scale  turbulence  for  the  Ameriflux  CO2  flux  network.  In  prin¬ 
ciple,  hot-wire  and  cold-wire  anemometers  enable  one  to  make  measurements  of  small-scale 
turbulence,  but  such  anemometers  are  unreliable,  often  suffering  mechanical  break-down. 
Furthermore,  they  disturb  the  flow  around  them.  Therefore,  there  is  a  need  for  a  new  genera¬ 
tion  of  sonic  anemometers/thermometers  with  increased  spatial  resolution.  The  main  goal  of 
the  present  paper  is  to  study  the  feasibility  of  using  a  sonic  anemometer  as  a  small  acoustic 
tomography  array  to  increase  the  spatial  resolution  in  measurements  of  temperature  and 
velocity. 

Ultrasound  travel-time  tomography  has  been  used  for  many  years  in  medicine  and  tech¬ 
nology  for  reconstructing  the  image  of  a  medium  through  which  the  ultrasound  waves  prop¬ 
agate.  In  these  fields,  the  number  of  ultrasound  propagation  paths  usually  greatly  exceeds 
the  number  of  unknown  parameters  of  the  medium  so  that  tomographic  reconstruction  is  an 
overdetermined  inverse  problem. 

The  first  experimental  implementation  of  acoustic  travel-time  tomography  of  the 
atmospheric  surface  layer  was  reported  by  Wilson  and  Thomson  (1994).  In  that  case,  the 
number  of  sound  propagation  paths  is  usually  less  than  the  number  of  temperature  and  wind 
velocity  values  to  be  reconstructed  at  selected  points  within  the  tomographic  region,  so 
that  one  deals  with  an  underdetermined  inverse  problem.  To  solve  this  problem,  the  time- 
dependent  stochastic  inversion  (TDSI)  algorithm  was  recently  developed  and  successfully 
applied  to  the  reconstruction  of  the  temperature  T  and  wind  velocity  v  fields  in  acoustic 
tomography  experiments  (Vecherin  et  al.  2006,  2008).  The  current  status  of  acoustic  tomog¬ 
raphy  of  the  atmospheric  surface  layer  is  reviewed  in  Ostashev  et  al.  (2009). 

Here,  we  suggest  that  one  consider  a  sonic  anemometer  as  a  small  acoustic  tomography 
array  and  use  appropriate  inverse  algorithms  for  reconstruction  of  the  fine  structure  of  the  T 
and  v  fields  inside  the  sonic  anemometer.  A  particular  modification  of  the  sonic  anemometer 
is  studied  when  the  number  of  its  transducers  is  doubled  and  the  TDSI  algorithm  is  used  for 
tomographic  reconstruction.  Numerical  simulations  of  temperature  and  velocity  measure¬ 
ments  with  the  sonic  anemometer  and  its  suggested  modification  are  implemented.  In  the 
simulations,  the  T  and  v  fields  inside  the  sonic  anemometer  and  its  modification  are  modelled 
with  quasi- wavelets.  The  quasi- wavelets  are  eddy-like  structures  of  temperature  and  velocity 
that  are  a  convenient  representation  of  turbulence  at  small  scales  (Wilson  et  al.  2009).  The 
numerical  simulations  allow  us  to  compare  the  errors  in  measurements  of  temperature  and 
velocity  with  the  sonic  anemometer  and  its  suggested  modification,  and  the  spatial  resolutions 
with  which  the  spectra  of  temperature  and  velocity  fluctuations  are  obtained. 
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The  paper  is  organized  as  follows.  In  Sect.  2,  the  geometry  and  principle  of  operation  of 
the  sonic  anemometer  and  its  suggested  modification  are  considered.  In  Sect.  3,  the  results  of 
numerical  simulations  are  presented  and  discussed.  Section  4  summarizes  the  results  obtained 
in  the  paper. 


2  The  Sonic  Anemometer  and  Its  Suggested  Modification 

2. 1  Sonic  Anemometer 

Although  the  design  of  sonic  anemometers  varies  substantially,  their  principle  of  operation 
is  similar  and  can  be  understood  by  considering,  for  example,  the  CSAT3  sonic  anemometer 
shown  in  Fig.  1.  There  are  three  transducers  at  both  the  “upper”  and  “lower”  levels  of 
the  instrument  that  point  downwards  and  upwards,  respectively.  In  numerical  simulations 
of  the  sonic  anemometer,  we  assume  that  the  upper  and  lower  transducers  are  separated 
in  the  vertical  direction  by  0.2  m.  These  transducers  form  three  sound  propagation  paths, 
marked  as  a,  b,  and  c,  whose  elevation  angles  0  with  respect  to  a  horizontal  plane  are  60°. 
The  projections  of  these  paths  on  the  horizontal  plane  located  in  the  middle  of  the  sonic 
anemometer  intersect  at  120°.  The  transducers  serve  alternately  as  sources  (transmitters)  and 
receivers.  As  a  result,  the  sound  pulses  propagate  in  the  opposite  directions  (i.e.,  upward 
and  downward)  of  each  path,  thus  enabling  measurements  of  six  travel  times  r ;  of  sound 
propagation  between  transducers.  Here,  i  =  1,  2, . . . ,  N,  where  N  =  6. 

A  sonic  anemometer  enables  measurements  of  the  path-averaged  values  of  the  sound 
speed  c  =  ^Jya  Ra  T  and  the  velocity  v  =  (vx,  vy,  vz).  Here,  ya  =  1.40  is  the  ratio  of 
specific  heats  for  dry  air,  Ra  =  287.058  J  kg-1  K-1  is  the  specific  gas  constant  for  dry 
air,  T  is  the  acoustic  virtual  temperature,  which,  hereafter,  is  referred  to  as  temperature,  and 
vx,  vy,  vz  are  the  velocity  components  in  a  Cartesian  coordinate  system  (x,y,z).  The  sound 
speed  c  and,  hence,  the  temperature  T  are  determined  from  a  sum  of  the  travel  times  T;  of 
sound  propagation  along  the  opposite  directions  of,  e.g.,  the  path  a.  The  projections  of  the 


Fig.  1  Schematic  of  the  CSAT3 
sonic  anemometer.  The 
ultrasound  propagation  paths 
a,  b,  and  c  are  shown  as  dashed 
lines  (reproduced  from  Horst  and 
Oncley  2006) 
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velocity  v  on  the  paths  a,  b,  and  c  are  determined  from  the  differences  of  the  travel  times  r ; 
of  sound  propagation  along  the  opposite  directions  of  these  paths,  respectively.  From  these 
projections,  it  is  straightforward  to  calculate  vx,  vy,  and  vz.  The  corresponding  relations  can 
be  found  in  Horst  and  Oncley  (2006),  and  in  the  CS  AT3  sonic  anemometer  instruction  manual 
(1998).  The  spatial  resolution  of  a  sonic  anemometer  is  determined  by  the  length  of  the  paths 
a,  b ,  and  c  shown  in  Fig.  1. 

2.2  A  Sonic  Anemometer  as  a  Tomography  Array 

A  suggested  modification  of  the  sonic  anemometer  is  shown  in  Fig.  2.  There  are  six  trans¬ 
ducers  at  both  the  upper  and  lower  levels  so  that  the  number  of  transducers  is  doubled  in 
comparison  with  a  standard  sonic  anemometer.  The  lower-  and  upper-level  transducers  are 
located  along  the  perimeter  of  a  square  with  a  side  length  of  0.1  m  and  are  separated  in  the 
vertical  direction  by  0.2  m,  i.e.,  the  same  distance  as  in  the  sonic  anemometer  considered 
above. 

In  Fig.  2,  each  transducer  works  alternately  as  a  source  and  receiver.  We  also  assume  that 
each  transducer  at  the  lower  level  can  receive  signals  from  every  transducer  at  the  upper 
level,  and  vice  versa.  This  increases  the  number  of  sound  propagation  paths  to  6  x  6  =  36. 
Taking  into  account  that  sound  signals  can  propagate  in  opposite  directions  of  each  path, 
the  total  number  of  travel  times  T;  that  can  be  measured  with  this  modification  of  the  sonic 
anemometer  is  N  =  72. 

2.3  Time-Dependent  Stochastic  Inversion 

In  the  considered  modification  of  the  sonic  anemometer,  reconstruction  of  the  temperature 
T (R,  t)  and  velocity  v(R,  t )  fields  is  done  with  the  TDSI  algorithm,  where  R  =  (x,y,z) 
are  the  Cartesian  coordinates  and  t  is  time.  It  is  convenient  to  express  these  fields  in  the 
following  form:  T( R,  t )  =  To(t)  +  7\(R,  t )  and  v(R,  t )  =  vq(0  +  vi(R,  t ),  where  the 
subscripts  “0”  and  “1”  correspond  to  the  mean  fields  (i.e.,  the  fields  averaged  over  the  tomo- 


Fig.  2  Schematic  of  the 
suggested  modification  of  the 
sonic  anemometer.  Transducers 
at  the  upper  (z  =  0.1  m)  and 
lower  (z  =  —0. 1  m)  levels  work 
alternately  as  ultrasound  sources 
and  receivers.  Green  lines  are  the 
ultrasound  propagation  paths 
between  transducers 
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graphic  volume)  and  their  fluctuating  parts,  respectively.  The  mean  fields  To(t)  and  vo(0 
are  reconstructed  from  the  measured  travel  times  r,  (t)  using  the  least-square  estimation. 
The  temperature  Tj  (R,  t)  and  velocity  vi  (R,  t)  fluctuations  are  reconstructed  with  the  TDSI 
algorithm. 

The  main  idea  of  the  TDSI  algorithm  is  to  assume  that  the  temperature  and  velocity  fluctu¬ 
ations  are  random  fields  with  the  known  spatial-temporal  correlation  functions  Bj  (R,  t)  and 
Bnm( R,  t),  respectively,  where  n,m  =  1,  2,  3.  The  cross-correlation  between  temperature 
and  velocity  fluctuations  is  assumed  to  be  zero,  which  is  justified  for  small-scale  turbulence. 
In  TDSI,  the  travel  times  T[  (tj)  are  repeatedly  measured  at  the  time  moments  tj ,  where 
j  =  1,  2, . . . ,  M,  and  then  used  in  the  reconstruction  of  the  temperature  and  velocity  fluctu¬ 
ations  at  a  specified  time. 

The  input  parameters  in  the  least-square  estimation  of  the  mean  fields,  Tq  and  Vo,  and  the 
TDSI  algorithm  for  reconstruction  of  the  fluctuations,  T\  and  Vi ,  are  the  travel  times  T;  (tj)  of 
sound  propagation,  the  distances  between  transducers,  and  parameters  of  the  spatial-temporal 
correlation  functions  Bj( R,  t)  and  Bnm( R,  t)  (the  latter  two  are  needed  for  TDSI).  The  total 
fields  T (R,  t)  =  To(t )  +  7\(R,  t )  and  v(R,  t )  =  vo(0  +  vi(R,  t),  and  the  estimated  errors 
in  their  reconstruction,  are  calculated  with  equations  derived  in  Vecherin  et  al.  (2006). 

Many  numerical  simulations  of  acoustic  tomography  of  the  atmosphere  were  implemented 
to  test  the  TDSI  algorithm.  In  the  simulations,  the  temperature  and  velocity  fields  within  a 
tomographic  area  or  volume  were  modelled  with  large-eddy  simulation  (LES),  e.g.,  Sullivan 
et  al.  (1994).  Different  forms  of  the  correlation  functions  of  temperature  and  velocity  fluc¬ 
tuations  were  used:  Gaussian,  exponential,  and  the  correlation  functions,  which  correspond 
to  the  von  Karman  spectra.  The  results  showed  that  the  accuracy  of  reconstruction  depends 
insignificantly  on  the  form  of  the  correlation  functions.  Therefore,  the  Gaussian  correlation 
functions  are  used  herein,  which  enable  some  integrals  in  the  TDSI  inversion  to  be  calculated 
analytically. 

The  Gaussian  spatial-temporal  correlation  function  of  statistically  homogeneous  and  sta¬ 
tionary  temperature  fluctuations  is  given  by  (Vecherin  et  al.  2006) 


Bj( R,  t )  =  Gj  exp 


(R  -  vo  tf 


(1) 


where  <7j  is  the  variance  of  temperature  fluctuations  and  Lj  is  the  characteristic  length. 
Equation  1  is  obtained  using  Taylor’s  frozen-turbulence  hypothesis.  The  spatial-temporal 
longitudinal  correlation  function  of  velocity  fluctuations,  Brr( R,  t ),  is  given  with  this  same 
equation  if  Gj  and  Lj  are  replaced  with  the  variance  g %  and  the  characteristic  length  Lv  of 
velocity  fluctuations,  respectively.  For  homogeneous  and  isotropic  turbulence,  the  spatial- 
temporal  correlation  functions  Bnm( R,  t )  are  expressed  in  terms  of  Brr( R,  t )  with  equa¬ 
tions  known  in  the  literature  (e.g.,  Eq.  6.34  in  Ostashev  1997).  The  characteristic  lengths 
Lj  and  Lv  determine  to  an  order  of  magnitude  the  smallest  scales  of  temperature  and 
velocity  fluctuations  that  can  be  resolved  in  reconstruction.  The  larger  scales  of  turbu¬ 
lence  reconstructed  with  TDSI  can  be  anisotropic  and  have  correlation  functions  that  dif¬ 
fer  from  the  Gaussian  correlation  functions.  For  example,  numerical  simulations  showed 
that  the  TDSI  algorithm  enables  accurate  reconstruction  of  the  LES  fields  corresponding 
to  stable,  neutral,  and  unstable  stratifications  (Ostashev  et  al.  2009).  There  is  a  trade-off 
between  the  values  of  Lj  and  Lv  and  the  expected  errors  in  the  reconstruction  of  tem¬ 
perature  and  velocity  fluctuations;  the  smaller  are  the  Lj  and  Lv  values  the  larger  are  the 
errors. 
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3  Numerical  Simulations 

3.1  Representation  of  Turbulence  as  Quasi- Wavelets 

In  numerical  simulation,  the  temperature  and  velocity  fields  have  been  modelled  as  quasi¬ 
wavelets  (Wilson  et  al.  2009).  Since  quasi- wavelets  are  not  yet  widely  used  in  the  literature, 
their  properties  are  briefly  reviewed  below. 

An  individual  temperature  quasi-wavelet  describes  a  fluctuation  in  temperature  with 
respect  to  the  mean  temperature, 

7T"(R)  =  qanf  (|R  -  ban\ //„)  ,  (2) 


where  /  is  a  spatially  localized  parent  function  (e.g.,  a  Gaussian  or  exponential  function), 
which  is  the  same  for  all  quasi- wavelets,  ban  is  the  centre  of  the  quasi-wavelet,  la  is  its 
size,  and  qan  is  its  amplitude.  The  index  a  =  1,  2,  . . . ,  N  corresponds  to  the  size  of  a 
quasi- wavelet,  with  l\  being  the  largest  size  and  In  the  smallest;  the  index  n  =  1,  2, . . . ,  Na 
describes  different  quasi-wavelets  with  the  same  size  la .  The  sum  of  quasi- wavelets  gives  us 
the  field  of  the  temperature  fluctuations  T\  (R).  To  construct  an  ensemble  of  quasi- wavelets, 
it  is  assumed  that  the  centres  ban  are  randomly  located  in  the  considered  volume,  the  mean 
value  of  the  amplitude  qan  is  zero,  and  the  second  moment  is  given  by  [qanq^m)  =  8ap8nmq^, 
where  8nm  is  the  Kronecker  delta  and  q %  is  the  mean  squared  amplitude.  The  remainder  of  the 
parameters  in  Eq.  2  are  deterministic  and  the  ratio  /a+i//a  =  e~ ^  is  assumed  constant.  The 
volume  occupied  by  the  quasi- wavelets  with  the  size  la  is  Va  =  Nal The  spatial  spectrum 
of  an  ensemble  of  quasi- wavelets  depends  on  the  form  of  a  parent  function,  parameter  /x, 
and  the  dependence  of  qa  and  Va  on  a.  With  these  choices,  realistic  spectra  of  temperature 
fluctuations  can  be  obtained,  such  as  the  Kolmogorov  and  von  Karman  spectra. 

Velocity  quasi-wavelets  are  constructed  similarly  (Wilson  et  al.  2009).  Thus,  quasi¬ 
wavelets  are  similar  to  wavelets  in  that  they  are  based  on  translations  and  dilations  of  a 
parent  function;  however,  their  positions  are  random.  Quasi-wavelets  enable  the  modelling 
of  intermittent,  inhomogeneous,  and  anisotropic  turbulence  with  prescribed  statistical  prop¬ 
erties.  One  of  the  applications  of  quasi- wavelets  is  to  generate  the  temperature  and  velocity 
fields  needed  in  numerical  simulations  of  sound  propagation  through  a  turbulent  atmosphere. 

In  this  paper,  quasi-wavelets  were  used  to  generate  statistically  homogeneous  and  isotropic 
temperature  and  velocity  fluctuations  inside  the  sonic  anemometer  volume  and  its  suggested 
modification.  The  smallest  and  largest  quasi-wavelets  were  chosen  as  l\  =  0.01  m  and  In  =  2 
m,  respectively,  while  other  parameters  of  the  quasi-wavelets  were  chosen  to  ensure  realistic 
power  spectra  of  temperature  and  velocity  fluctuations. 

In  numerical  simulation,  the  size  of  the  computational  grid  along  the  x,  y,  and  z-axes  is 
equal  to  the  smallest  quasi-wavelet,  i.e.,  0.01  m.  The  quasi-wavelet  temperature  and  veloc¬ 
ity  fields  were  created  in  a  volume  of  20m  x  0.1  m  x  0.2m;  in  the  simulations,  this  volume 
moved  in  the  direction  of  the  v-axis  with  velocity  U  =  2 ms-1 ;  the  total  time  of  the  numer¬ 
ical  experiment  was  (20m)/(2ms_1)  =  10  s.  Figure  3  shows  a  snapshot  of  the  temperature 
and  three  components  of  the  velocity  fields  inside  the  sonic  anemometer  and  its  suggested 
modification  obtained  with  quasi- wavelets.  In  the  numerical  simulations,  the  discretization 
time  was  0.005  s  so  that  the  total  number  of  snapshots  of  the  T  and  v  fields  was  2000.  For 
each  snapshot,  the  travel  times  r*  along  all  sound  propagation  paths  in  the  sonic  anemometer 
(Fig.  1)  and  its  suggested  modification  (Fig.  2)  were  calculated. 
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Fig.  3  Snapshot  of  the  temperature  and  three  components  of  the  velocity  fields  obtained  with  quasi- wavelets. 
Upper  left :  T(x,y,z),  upper  right :  vx(x,  y,  z),  lower  left :  vy(x,  y,  z),  and  lower  right :  vz(x,  y,  z) 
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3.2  Reconstruction  of  Temperature  and  Velocity  Fields 

The  travel  times  r i  for  the  sonic  anemometer  were,  then,  used  to  calculate  the  path- 
averaged  values  of  temperature  and  velocity  with  equations  known  in  the  literature 
(e.g.,  Kaimal  and  Finnigan  1994).  The  obtained  time  series  of  T ,  vx,  vy,  and  vz  are  depicted  as 
red  lines  in  the  left-hand  plots  in  Figs.  4,  5,  6  and  7,  respectively;  the  blue  lines  are  the  quasi¬ 
wavelet  values  of  T,  vx,  vy,  and  vz  at  the  centre  of  the  sonic  anemometer.  The  agreement 
between  these  quasi-wavelet  values  and  the  path-averaged  values  of  T,  vx,  vy,  and  vz  recon¬ 
structed  in  numerical  simulations  of  the  sonic  anemometer  is  relatively  good.  This  comparison 
shows  to  what  extent  the  path-averaged  values  of  temperature  and  velocity  represent  the  cor¬ 
responding  values  at  the  centre  of  a  sonic  anemometer.  The  root-mean- square  errors  ( RMSEs ) 
between  the  red  and  blue  lines  in  the  left  plots  in  Figs.  4-7  are  Oj  =  0.09  °C,  or®  =  0.086 
ms-1,  =  0.083  ms-1,  and  =  0.056  ms-1,  respectively,  where  the  superscript 

“s”  stands  for  “sonic”.  The  RMSE  for  the  vertical  velocity,  or  ®  ,  is  about  1.5  times  smaller 
than  the  RMSE  for  the  horizontal  components  of  the  velocity,  or®  and  a ®  ,  and  is  due  to  the 
fact  that  the  sound  propagation  paths  in  the  sonic  anemometer  have  relatively  large  elevation 
angles  0  =  60°.  Indeed,  using  equations  known  in  the  literature,  e.g.,  Eqs.  9a-9c  in  Horst  and 
Oncley  (2006),  vx,  vy,  vz  are  calculated  from  the  velocity  components  va,  Vb,  vc  measured 
with  the  sonic  anemometer  along  the  sound  propagation  paths  a,  b,  c  shown  in  Fig.  1 .  If  Av 
is  the  error  in  the  measurement  of  va,  Vb,  vc,  it  follows  from  these  equations  that  the  error  in 
the  measurement  of  vx  or  vy  is  proportional  to  Av/  cos  0,  while  the  error  in  the  measurement 
of  vz  is  proportional  to  A  v/  sin  0.  For  0  =  60°,  the  former  error  is  larger  than  the  latter  by 
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Fig.  4  Time  series  of  temperature,  T  (t).  Blue  lines  are  the  quasi- wavelet  values  of  temperature  at  the  centre 
of  the  sonic  anemometer  and  its  suggested  modification.  Red  lines  are  the  values  of  temperature  reconstructed 
in  numerical  simulations  of  the  sonic  {left  plot )  and  its  modification  (right  plot ) 


Fig.  5  Time  series  of  the  x-component  of  the  horizontal  velocity,  vx(t).  Blue  lines  are  the  quasi- wavelet 
values  of  vx  at  the  centre  of  the  sonic  anemometer  and  its  suggested  modification.  Red  lines  are  the  values  of 
vx  reconstructed  in  numerical  simulations  of  the  sonic  ( left  plot )  and  its  modification  ( right  plot ) 


tan  60°  ~  1.73,  which  is  close  to  the  ratios  /a ^  =  1.54  and  a*  /rr^  =  1.48  obtained  in 

numerical  simulations. 

Then,  the  travel  times  r i  for  the  considered  modification  of  the  sonic  anemometer  were 
used  as  input  parameters  for  reconstruction  of  the  mean  and  fluctuating  parts  of  the  tem¬ 
perature  and  velocity  fields.  In  the  TDSI  algorithm  for  reconstruction  of  the  fluctuations, 
three  consecutive  values  of  the  travel  times  t;  were  employed.  Furthermore,  the  square  roots 
of  the  variances  of  the  spatial-temporal  correlation  functions  of  temperature  and  the  three 
components  of  velocity  fluctuations  were  calculated  for  the  considered  quasi-wavelet  field, 
and  were  0.074°C,  0.057  ms-1,  0.058  ms-1,  and  0.064  ms-1,  respectively;  the  charac¬ 
teristic  length  of  these  fluctuations  was  Lp  =  Lv  =  0.03  m.  In  the  right-hand  plots  in 
Figs.  4-7,  the  time  series  of  the  reconstructed  values  of  T ,  vx,  vy,  and  vz  in  the  centre  of  the 
considered  modification  of  the  sonic  anemometer  are  shown  as  red  lines,  respectively.  The 
blue  lines  are  again  the  quasi-wavelet  values  of  T,  vx,  vy ,  and  vz  in  the  centre  of  the  sonic 
anemometer,  i.e.,  the  blue  lines  in  the  left-hand  and  right-hand  plots  in  these  figures  are  the 
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Fig.  6  Time  series  of  the  y-component  of  the  horizontal  velocity,  vy(t).  Blue  lines  are  the  quasi- wavelet 
values  of  vy  at  the  centre  of  the  sonic  anemometer  and  its  suggested  modification.  Red  lines  are  the  values  of 
vy  reconstructed  in  numerical  simulations  of  the  sonic  ( left  plot)  and  its  modification  ( right  plot) 
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Fig.  7  Time  series  of  the  vertical  component  of  the  velocity,  vz(t).  Blue  lines  are  the  quasi-wavelet  values 
of  vz  at  the  centre  of  the  sonic  anemometer  and  its  suggested  modification.  Red  lines  are  the  values  of  vz 
reconstructed  in  numerical  simulations  of  the  sonic  ( left  plot)  and  its  modification  ( right  plot) 


same.  It  follows  from  the  right-hand  plots  in  Figs.  4-7  that  the  agreement  between  the  quasi¬ 
wavelet  and  reconstructed  values  of  T,  vx,  vy,  and  vz  is  very  good,  and,  in  fact,  is  better  than 
for  the  sonic  anemometer  (see  left-hand  plots).  The  RMSEs  between  the  quasi- wavelet  and 
reconstructed  values  of  T,  vz,  vx ,  and  vy  are  smaller  than  those  for  the  sonic  anemometer: 
a™  =  0.75 a*  =  0.067  °C,  a,m  =  0.91a*  =  0.077  m  s“\  a,"1  =  0.85a*  =  0.07  m  s“\ 

1  1  UX  UX  7  Vy  Vy 

and  or™  =  0.62(7^  =  0.035  m  s-1,  respectively.  Here,  the  superscript  “m”  stands  for  the 
modification  of  the  sonic  anemometer. 

The  results  of  the  numerical  simulations  shown  in  Figs.  4-7  clearly  indicate  that  the  sug¬ 
gested  modification  of  the  sonic  anemometer  enables  a  more  accurate  reconstruction  of  the 
temperature  and  vertical  velocity,  and  somewhat  better  reconstruction  of  the  horizontal  veloc¬ 
ity.  These  results  are  used  in  the  next  sub-section  to  study  the  reconstruction  of  the  spectra  of 
temperature  and  velocity  fluctuations  with  the  sonic  anemometer  and  its  modification.  Note 
that  the  TDSI  algorithm  was  also  applied  to  the  reconstruction  of  temperature  and  velocity 
fields  inside  the  sonic  anemometer  (Fig.  1).  Though  not  presented  here  for  brevity,  the  results 
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obtained  showed  no  significant  improvement  in  the  reconstruction  of  temperature  and  the 
three  components  of  the  velocity  in  comparison  with  their  path- averaged  values. 

3.3  Spectra  of  Temperature  and  Velocity  Fluctuations 

The  time  series  T,  vx,  vy,  and  vz  obtained  in  numerical  simulations  of  the  sonic  anemometer 
and  its  suggested  modification  and  the  corresponding  values  of  the  quasi-wavelets  in  the 
centre  of  the  sonic  were  employed  to  calculate  the  mean  values  of  temperature  and  three 
components  of  the  velocity  and  their  fluctuating  parts.  The  fluctuations  were  then  used  to 
calculate  the  power  spectra  SSj(k),  S™  (k),  and  Sjw  (k),  where  “s”  stands  for  the  sonic,  “m” 
for  its  suggested  modification,  “qw”  for  quasi- wavelets,  and  k  is  the  turbulence  wavenumber. 
The  subscript  j  stands  for  T,  x,  y,  or  z  and  corresponds  to  fluctuations  in  temperature  or  the 
three  components  of  the  velocity. 

The  spectral  transfer  functions  are  defined  as  the  ratios  of  the  power  spectra  (e.g.,  Horst 
and  Oncley  2006;  Angelou  et  al.  2012),  viz. 


(3a) 


(3b) 


The  spectral  transfer  functions  characterize  the  effects  of  path  averaging  and  errors  in  recon¬ 
struction  of  the  power  spectra  of  temperature  and  velocity  fluctuations  obtained  in  numerical 
simulations  of  the  sonic  anemometer  and  its  suggested  modification.  The  spectral  transfer 
functions  are  plotted  in  Fig.  8;  the  red  and  blue  lines  correspond  to  Hj  (k)  and  respec¬ 

tively.  For  a  perfect  reconstruction  of  the  spectra,  these  lines  would  coincide  with  the  black 
horizontal  lines  corresponding  to  the  case  Hj(k)  =  =  1. 

It  follows  from  the  upper  left-hand  and  lower  right-hand  plots  in  Fig.  8  that  the  sonic 
anemometer  enables  the  correct  reconstruction  of  both  the  temperature  and  vertical  velocity 
spectra  for  the  wavenumbers  k  smaller  than  about  k^s  =  10  m-1,  where  “ps”  stands  for 
the  power  spectrum.  Let  us  compare  k^s  with  a  similar  value  obtained  by  Horst  and  Oncley 
(2006)  who  studied  analytically  the  effect  of  path  averaging  in  the  CS  AT3  sonic  anemometer. 
It  follows  from  Fig.  3  in  that  paper  that  the  spectral  transfer  functions  Hj(k)  are  close  to  1 
for  k  smaller  than  about  k g  =  1/L,  where  L  is  the  distance  between  the  transducers.  In  the 
sonic  anemometer  considered  herein,  L  =  0.23  m  so  that  k q  =  4.3  m_1.  The  difference 
between  k q  and  k^s  can  be  explained  by  the  variability  in  the  red  (and  blue)  curves  in  Fig.  8, 
which  is  due  to  the  errors  in  reconstruction  and  which  makes  it  difficult  to  determine  the 
exact  value  of  the  cut-off  frequency  k^s.  It  also  follows  from  the  upper  left-hand  and  lower 
right-hand  plots  in  Fig.  8  that  the  suggested  modification  of  the  sonic  anemometer  enables 
the  correct  reconstruction  of  the  temperature  and  vertical  velocity  spectra  for  k  smaller  than 
about  k™  =  100  rad  m-1,  thus  increasing  its  spatial  resolution  by  a  factor  of  k™/k^s  =  10. 

The  spectral  transfer  functions  for  the  horizontal  velocity  components,  vx  and  vy,  recon¬ 
structed  in  numerical  simulations  of  the  sonic  anemometer  and  its  modification  are  shown  in 
the  upper  right-hand  and  lower  left-hand  plots  in  Fig.  8,  respectively.  It  follows  from  these 
figures  that  the  spectral  transfer  functions  significantly  deviate  from  1  in  the  range  k  >  10 
m-1  and  that  the  reconstruction  accuracy  is  about  the  same  for  the  sonic  anemometer  and 
its  modification.  These  results  are  due  to  predominantly  vertical  orientations  of  the  sound 
propagation  paths  for  both  the  sonic  anemometer  and  its  modification,  when  the  horizontal 
components  of  the  velocity  are  reconstructed  less  accurately  than  the  vertical  component. 
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Fig.  8  Spectral  transfer  functions  defined  through  the  power  spectra  and  obtained  in  numerical  simulations 
of  the  sonic  anemometer  (Hj(k),  red  lines )  and  its  suggested  modification  (H™(k),  blue  lines).  The  black, 
horizontal  lines  correspond  to  the  case  Hj(k)  =  H™  (k)  =  1.  Upper  left  T  fluctuations;  upper  right  vx 
fluctuations;  lower  left  vy  fluctuations;  and  lower  right  vz  fluctuations 


Figure  8  also  shows  that  the  fastest  decrease  of  the  spectral  transfer  functions  obtained 
in  numerical  simulations  of  the  sonic  anemometer  is  for  the  temperature  power  spectrum, 
followed  by  the  power  spectra  of  vz,  vy,  and  vx.  The  spectral  transfer  functions  plotted  in 
Fig.  3  in  Horst  and  Oncley  (2006)  have  qualitatively  the  same  behaviour. 

Angelou  et  al.  (2012)  argue  that  the  effect  of  path  averaging  on  measurements  of  the 
radial  component  of  the  wind  velocity  with  a  laser  anemometer  should  be  taken  into  account 
using  cross-spectra  rather  than  power  spectra.  Indeed,  the  noise  in  the  experimental  data 
can  be  amplified  in  the  power  spectrum,  while  it  is  cancelled  in  the  cross-spectrum.  In 
the  present  paper,  the  noise  in  the  data  is  due  to  the  errors  in  numerical  reconstruction  of 
the  temperature  and  velocity  fields.  As  a  result,  the  spectral  transfer  functions  in  Fig.  8 
have  noticeable  variability  and  are  >1  for  some  values  of  k  though  they  should  always 
be  <1. 

Following  Angelou  et  al.  (2012),  the  cross-spectra  xj’qw(k)  between  fluctuations  in  tem¬ 
perature  or  velocity  components  obtained  in  numerical  simulations  of  the  sonic  anemometer 
and  the  corresponding  values  of  the  quasi-wavelets  in  the  centre  of  the  sonic  were  calcu¬ 
lated.  Similarly,  the  cross-spectra  ^>qw(£)  between  fluctuations  in  temperature  or  velocity 
components  obtained  in  numerical  simulations  of  the  suggested  modification  of  the  sonic 
anemometer  and  the  values  of  the  quasi- wavelets  were  also  calculated.  The  spectral  transfer 
functions  are  now  defined  similarly  to  those  in  Angelou  et  al.  (2012), 
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Fig.  9  Spectral  transfer  functions  defined  through  the  cross-spectra  and  obtained  in  numerical  simulations 
of  the  sonic  anemometer  ( hSj(k ),  red  lines )  and  its  suggested  modification  (h^ik),  blue  lines).  The  black, 
horizontal  lines  correspond  to  the  case  hSj(k)  =  h^ik)  =  1.  Upper  left  T  fluctuations;  upper  right  vx 
fluctuations;  lower  left  vy  fluctuations;  and  lower  right  vz  fluctuations 


The  spectral  transfer  functions  hs-(k)  and  h™(k)  are  plotted  in  Fig.  9;  the  red  and  blue 
lines  correspond  to  hSj(k)  and  h™(k),  respectively;  the  black  horizontal  lines  to  the  case 
h*j(k)  =  hj(k)  =  1. 

The  spectral  transfer  functions  in  Fig.  9  confirm  conclusions  made  from  Fig.  8:  the  sug¬ 
gested  modification  of  the  sonic  anemometer  enables  us  to  increase  the  spatial  resolution  of 
the  temperature  and  vertical  velocity,  while  the  spatial  resolution  of  the  horizontal  velocity 
remains  approximately  unchanged.  These  results  are  actually  more  evident  in  Fig.  9  since  the 
cross-spectra  in  the  upper  left-hand  and  lower  right-hand  plots  in  Fig.  9  have  less  variability 
than  similar  power  spectra  in  Fig.  8.  Furthermore,  the  spectral  transfer  functions  in  Fig.  9  are 
^  1,  as  they  should  be.  It  follows  from  the  upper  left-hand  and  lower  right-hand  plots  in  this 
figure  that  the  sonic  anemometer  enables  the  correct  reconstruction  of  both  the  temperature 
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and  vertical  velocity  spectra  for  wavenumber  k  smaller  than  about  k^s  =  5  m_1,  where  the 
subscript  “cs”  stands  for  cross-spectra.  This  value  of  /^s  is  close  to  =  4.3  m-1  obtained 
from  the  analytical  analysis.  Furthermore,  it  follows  from  these  plots  that  the  spectral  trans¬ 
fer  functions  for  temperature  and  vertical  velocity  obtained  in  numerical  simulations  of  the 
suggested  modification  of  the  sonic  anemometer  are  close  to  1  for  k  <  50  m-1.  Thus,  the 
sonic  anemometer  as  a  small  tomography  array  enables  us  to  increase  the  spatial  resolution 
by  a  factor  of  about  50/5  =  10.  For  k  >  80  m-1,  the  spectral  transfer  functions  in  Fig.  9 
increase  or  oscillate  due  to  numerical  instability  caused  by  small  values  of  the  power  spectra 
and  cross-spectra  at  such  wavenumbers. 


4  Conclusions 

To  increase  the  spatial  resolution  of  a  sonic  anemometer,  it  is  suggested  that  one  consider  it  as 
a  small  acoustic  tomography  array  and  apply  appropriate  inverse  methods  for  reconstruction 
of  the  temperature  and  velocity  fields.  A  particular  modification  of  the  sonic  anemometer  was 
studied  when  the  number  of  its  transducers  is  doubled  and  reconstruction  of  the  temperature 
and  velocity  fields  is  done  with  the  TDSI  algorithm.  Numerical  simulations  of  the  sonic 
anemometer  and  its  suggested  modification  were  implemented  with  the  temperature  and 
velocity  fields  modelled  as  quasi- wavelets.  The  simulations  showed  that  the  sonic  anemome¬ 
ter  enables  measurements  of  the  time  series  of  temperature  and  velocity  with  a  relatively  small 
root-mean- square  error  ( RMSE ).  The  suggested  modification  of  the  sonic  anemometer  allows 
even  more  accurate  measurements,  especially  for  temperature  and  the  vertical  component  of 
velocity,  when  the  corresponding  RMSE  is  smaller  than  that  for  the  sonic  anemometer  by 
factors  0.75  and  0.62,  respectively.  The  increase  in  accuracy  is  less  pronounced  for  the  hori¬ 
zontal  components  of  velocity  due  to  the  predominantly  vertical  orientations  of  the  ultrasound 
propagation  paths  in  the  sonic  anemometer  and  its  modification. 

Another  advantage  of  the  considered  modification  of  the  sonic  anemometer  is  a  signif¬ 
icant  increase  in  the  spatial  resolution  of  the  spectra  of  temperature  and  vertical  velocity 
fluctuations.  The  numerical  simulations  showed  that  the  sonic  anemometer  enables  correct 
measurements  of  these  spectra  up  to  wavenumbers,  that  is  about  the  inverse  of  the  distance 
between  its  transducers.  This  result  clearly  indicates  the  effect  of  path  averaging  on  the  mea¬ 
surements  of  the  spectra  and  agrees  with  that  obtained  by  Horst  and  Oncley  (2006)  who 
used  an  analytical  approach.  On  the  other  hand,  the  considered  modification  of  the  sonic 
anemometer  enables  correct  reconstruction  of  these  spectra  up  to  the  wavenumbers  that  are 
about  ten  time  larger  than  those  for  the  sonic  anemometer.  Its  spatial  resolution,  generally, 
depends  on  the  inverse  algorithm  for  reconstruction  of  the  temperature  and  velocity  fields 
and  also  improves  by  increasing  the  number  of  transducers. 

From  a  technological  point  of  view,  the  considered  modification  of  the  sonic  anemometer 
seems  feasible;  other  modifications  might  also  be  feasible  with  more  transducers  at  the  upper 
and  lower  levels.  With  the  increased  number  of  transducers,  the  spatial  resolution  can  be 
increased  further.  Also,  the  inverse  problem  in  the  reconstruction  of  temperature  and  velocity 
might  become  overdetermined  so  that  standard  inverse  algorithms  (e.g.,  the  least-square 
estimation)  can  be  employed.  Note  that  the  increased  number  of  transducers  might  increase 
the  flow  distortion.  However,  the  distance  between  transducers  can  be  increased  so  that  the 
temperature  and  velocity  fields  in  the  centre  of  a  sonic  volume  will  not  be  distorted  and  can 
be  reconstructed  with  the  use  of  acoustic  tomography.  Since  acoustic  tomography  provides 
turbulent  fields  at  all  points  inside  a  sonic  anemometer,  it  can  also  be  used  to  study  the  flow 
distortion  caused  by  transducers. 
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